Phase separation in the combined Falicov-Kimball and static Holstein model 
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The Falicov-Kimball model has been known to phase separate when the correlation strength is 
large and the temperature is low. We examine the stability of phase separation under the influence of 
an additional electron-phonon interaction that tries to force the electron density to be homogeneous. 
We find that the electron-phonon interaction eventually destroys the phase-separated state once it 
becomes larger than some critical value. The results of this work may have an influence on theories 
of collossal magentoresistance materials. 
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The phenomenon of strong electron correlations driv- 
ing phase separation has been long conjectured in the 
Falicov-Kimball modeltl and in the Hubbard modelo. 
Little progress has been made on rigorously establish- 
ing phase separation in the Hubbard model, but it has 
been proved to exist in the Falicov-Kimball model in 
one dimensiorio, in infinite dimensions!], and in all other 
dimensions™]. In spite of these rigorous results, there has 
still been only limited calculations of phase diagrams and 
the transition temperature as a function of the parame- 
ters in the model. Essentially all finite-temperature cal- 
culationSphave been performed in the infinite-dimensional 
limiijj'LlQQ. Here we continue with this tradition, and il- 
lustrate how competing (finite) interactions can destroy 
the phase-separated state, even if the Falicov-Kimball in- 
teraction strength is infinitely large. 

The system we choose to study is the harmonic 
versioit,|Of the spinless Falicov-Kimball-static-Holstein 
.defie. 



model 



The Hamiltonian is 



(id) 



9 K E^ + E^ ~ V)(4°i "As), 



(1) 



where cj (cj) is the creation (annihilation) operator for an 
itinerant electron at site i, fj (fa) is the creation (annihi- 
lation) operator for a localized electron at site i, t* j1\fd is 
the hopping integral between nearest-neighbor sites (each 
pair is summed over twice for hermiticity) , Ef is the local 
f-clcctron site energy, U is the Falicov-Kimball interac- 
tion, Xi is the static phonon coordinate, k is the spring 
constant (chosen to be 1 here), and g is the electron- 
phonon interaction. The symbol p e denotes the aver- 
age itinerant-electron concentration (ere}. We examine 
the system in the infinite-dimensional limiui3 and take 
t* as our energy unit. A common chemical potential /i 
is used to fix the total electron concentration (itinerant 
plus localized). Note that we have chosen the phonons 
to interact with the local charge fluctuation, which is a 
common choice, but not the only one that could be made. 



Since the term that is subtracted depends on the local- 
ized electron concentration it is not just a constant that 
can be included or removed; it has a nontrivial effect on 
the value of the Hclmholz free energy, and thereby on the 
phase separation. For concreteness we solve our system of 
equations on an infinite-coordination Bcthe lattice, which 
has a semicircular electronic density of states 



Z7T 



(2) 



The system is solved by employing dynamical mean- 
field theory. The procedure simply combines the strategy 
used for the Falicov-Kimba]Lmodclt2l with that used for 
the static Holstein modelE-alij. We summarize the pro- 
cedure here to standardize our notation. Our interest is 
in the U — > oo limit. In that case, the itinerant elec- 
trons and the localized electrons are not allowed to sit 
at the same lattice site. The dynamical mean-field the- 
ory is straightforward to evaluate. The partition function 
becomes 



Z = 



dx2e-I^Y[ G ° 1{iuJn) - 9X 



ioj n 



1 



x exp{-/3[-(gx - p)p e + -kx ]} 

' 00 dx2e~^-^+u) TT G^(iu n )-gx-U 

x exp{-/3{~(gx- p)p e + E f - p + -kx 2 }}, (3) 

with P — l/T the inverse temperature, Go(iu) n ) = 
l/(iu) n + p. — X n ) the effective medium, iuj n = ii:T(2n+l) 
the Fermionic Matsubara frequency, and \(ioj n ) = X n the 
dynamical mean field. Once the partition function is de- 
termined (in the limit U — > oo), then one can calculate 
the single-particle electronic Green's function via 



G(iu> n ) = G n 



dx- 



w (x) 



(4) 



G^iiuJn) - gx' 
with wq(x) the phonon distribution function (when no 



2 



localized-electron occupies the lattice site), which satis- 
fies 



wo(x) — 2e 



|(sx-m) TT G o 1 ( iu >n) - gx 



x exp{-/3[-( 5 a; - p)p e + -nx 2 ]}/Z; (5) 
the integral of wq(x) satisfies a sum rule 

/oo 
dxwo(x). (6) 
-oo 

[The sum rule is not one, because we have taken the limit 
U — * oo. If C7 is finite, we would have a second phonon 
distribution function W\(x) to consider when calculating 
the Green's function, which enters when the site is oc- 
cupied by a localized-electron. In the U — oo limit, we 
can always rewrite wq(x) = P(x) — w\{x) where the inte- 
gral of P(x) is one and can be interpreted as the U = oo 
phonon distribution function, and w\{x) integrates to the 
localized electron concentration. But in this work, we 
never distinguish these objects, so there is no need to 
elaborate further on them.] To complete the loop for the 
self-consistent solution of the problem, we need to note 
that Dyson's equation is 



and 



G,. 



f?o 1 (^„) = G- 1 +S r , 



<M e ) -■ — ; v ' 

_ 2 lU) n + p, - L„ - e 



(7) 



(8) 



with £(iw n ) = Ti n the itinerant-electron self energy. This 
provides a closed set of equations that is needed to solve 
the many-body problem. 

The computational algorithm that weJpllow is the it- 
erative procedure introduced by Jarrelllij: (i) we begin 
with the self energy set equal to zero; (ii) Eq. (||) is used 
to determine G n ; (iii) Eq. (0) is used to extract the ef- 
fective medium G ; (iv) Eqs. (^), (g), and (|4|), then de- 
termine G„; and (v) Eq. (0) is employed to extract the 
new self energy. Then steps (ii-v) are repeated until con- 
vergence is reached. Once the equations have converged 
on the imaginary axis, one can take the function wo(x) 
and the analytic continuation of Eqs. (|J), (|7|), and (|4|) to 
determine the Green's function on the real axis as well. 

In order to study phase separation, we also need to cal- 
culate the free energy. This is accomplished by following 
the procedure of Brandt and MielschQ. The Helmholz 
free energy satisfies 



-T lnZ + ^lnG„ 

n 

Y, f dep{e)\n . 1 
^ J - 2 iuj n + (i- E n - e 

(E f -p)p f . J (9) 



It is important to include a sufficient number of Matsub- 
ara frequencies in the summation to achieve convergence 
of the summation; indeed this ultimately limits the ac- 
curacy of our results. 

It is well known that if both p e and pf lie between 
zero and one, then the ground state is phase separatedu. 
The key point in this work is that we fix Ef. If the 
chemical potential lies below Ef then the ground state 
contains only itinerant electrons and is not phase sepa- 
rated. As the chemical potential is increased, we reach 
a critical value of p where the average localized elec- 
tron filling becomes nonzero, and the system goes into 
a phase-separated state. It is straightforward to deter- 
mine the value of this critical chemical potential when 
g = 0. The interacting density of states is just the semi- 
circle in Eq. (|J) when there are no localized electrons and 
it becomes 



p™t(e) = ^V 4 ! 1 - i'f> -<- 



(10) 



when there are pj interacting electrons. So we must 
compare the energy with no localized electrons E gs — 
— (4— /^ 2 ) 3//2 /67r to the mixture of states with all localized 
particles and no itinerant particles (with weight pf) and 
the state with no localized particles but with an electron 
filling of p e = ptot - (1 - Ptot)pf (with weight I — Pf) in 
the limit where pf approaches zero. Note that one must 
adjust the chemical potential p* of the mixed state, to 
have the correct electron density. After performing some 
straightforward algebra, we find a transcendental equa- 
tion for the critical value of the chemical potential 



E f = - 



67T 



1 . _i A 4 
— sin — 

TT 2 



(11) 



As Ef ranges from —2 to 0, the critical total electron 
density for phase separation ranges from to 1; for Et > 
the system does not phase separate as T — > 0. 

We will examine the case with Et = —0.5 here. 
Eq. @ is solved by p = -0.145 and p e + Pf = 0.454 
for the lower boundary where phase separation begins. 
Hence phase separation enters when p tot > 0.454. To 
determine the phase separation, we must perform a 
Maxwell construction. If we plot the electron filling as a 
function of the chemical potential, we will find a multival- 
ued function in the region of phase separation. One way 
to perform the Maxwell construction is by drawing a line 
that has equal areas above and below, with the points of 
intersection corresponding to the densities of the phases 
that compose the mixture of the phase-separated state 
(the other way is to construct the convex hull of the free- 
energy curve; both methods yield the same results). The 
double- valued function and the Maxwell construction are 
plotted for g = 1.2 and T = 0.04 in Fig. |. It is nontrivial 
to generate such a curve. The technical difficulty arises 
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FIG. 1: The chemical potential is plotted versus the total 
electron filling when g = 1.2, E f = -0.5 and T = 0.4. Note 
the double-valued chemical potential indicating phase sepa- 
ration. The Maxwell construction is also shown in the figure 
by the dashed line. The total area of the solid curve above 
and below the dashed line is equal to zero. The dotted lines 
indicate the total electron densities of the states that form 
the mixture. 



from the fact that for a fixed value of Ef, there are mul- 
tiple solutions of pf which have the same ptot when one is 
in the phase separated region. Hence, we sometimes need 
to either adjust pf and calculate Ef once the equations 
have converged, or employ a two-dimensional rootfinder 
for p and pf to get the equations to converge. Even so, 
we find regions of parameter space where we are unable 
to stabilize the algorithm. 

We examine the transition temperature for different 
electron-phonon coupling strengths in Fig. |^. The phase- 
separated region hardly changes at zero temperature as 
the electron-phonon coupling increases, but the transi- 
tion temperatures are sharply reduced. We run into nu- 
merical problems with the stability of our iterative so- 
lutions when g is increased to 1.2 and above. We are 
unable to find any strategy to solve for the phase sep- 
aration at low temperature in this region of parameter 
space. We expect, however, that the phase separation 
continues to be suppressed, and conjecture that there is a 
critical value of the electron-phonon coupling where the 
phase separation disappears, but we cannot determine 
it's value with our numerical algorithms. In addition 
to the binodal (first-order) phase transition temperature, 
we also include the spinodal decomposition temperature 
(below which the system cannot be supercooled) for the 
g = case. This is found by determining where the 
k = charge susceptibility diverges. Of course the bin- 
odal and spinodal temperatures must agree at the peak 
of the curve. 

The disappearance of the phase separation occurs 



when the electron-phonon interaction energy becomes 
similar in magnitude to the kinetic-energy gain of the 
phase separated state. This occurs at a finite value of 
g because the kinetic energy gain is finite, and does not 
diverge as U — ► oo. 

To understand this phenomenon further, we examine 
in detail the system in the low-temperature limit. We 
already saw that when g = there is a minimum total 
electron density that must be reached before phase sepa- 
ration occurs. Similarly, when g ^ 0, we find that such a 
critical density continues to hold for small g. This occurs 
because, as T is lowered, the function wq(x) becomes 
sharply peaked at x = 0, eventually becoming a delta 
function at T = 0. Hence, the system remains a Fermi 
gas in this regime of parameter space. As g is increased to 
a critical value, wq(x) develops a local minimum at x = 
and two maxima, one at x < and one at x > 0. At 
this point, the ground state becomes a non-Fermi-liquid, 
and we believe that the non-Fermi-liquid character sup- 
presses the phase separation, ultimately eliminating it. 
To be more quantitative, we compute 



n riln wq(x) 
dx 



9(Pe 



1 



G (iu> n ) - gx 



When this derivative vanishes, we have an extremum of 
the phonon distribution function. When pf = 0, this 
equation always has a solution at x = because the 
sum over Matsubara frequencies equals p e — 1/2. If g 
is increased, we reach a point where the second deriva- 
tive becomes positive, and the x = solution becomes 
a local minimum instead of a maximum. At this point 
two nonzero x solutions exist, and the system becomes a 
non-Fermi-liquid. The second derivative (at x — 0) is 



T 



d 2 In w (x = 0) 
dx 2 



g 



•■Tj2G 2 (iu; n ), 



(13) 



which becomes 

n d 2 \nwo(x = 0) 



T- 



dx 1 



It" J lhp\l ;I{('C"(( i 

,2\ " 



-k + g- 



12tt 



(14) 



in the low-temperature limit. For p = —0.145 the critical 
value of g to enter the non-Fermi-liquid state is 2.18. 

It is useful to discuss the impact of this work on 
models for the collossal magnetoresistance manganite 
materials. There are a number of different elements 
that appear to be important in describing the CMB, 
materials. All models start with double-exchangeO, 
but then need something else to enhance the metal- 
insulator transition. Three main ideas have emjerged 
to do this: a strong electatfwpkonon interactiontftE3E3; 
strong disorder scatteringE3E3 , Lil; and proximity to phase 
separations. The results of this investigation indicate 
that the three elements are not fully independent. In- 
deed, strong disorder scattering can lead to phase separa- 
tion, and strong electron-phonon scattering can suppress 
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FIG. 2: The first-order phase transition as a function of 
the total electron filling (the solid dots show the results of the 
Maxwell construction). We examine four cases g = 0, g = 0.6, 
g = 0.8, and g — 1.0. The dashed line is the spinodal decom- 
position temperature (found by calculating the divergence of 
the k — charge susceptibility) for the g — case. 



the phase separation. One of the consequences of this 
work is that one needs to carefully investigate the im- 
plications of all three of these physical effects, and how 
they interrelate, to construct a meaningful theory for the 
manganites. 

In conclusion, we have analyzed the suppression of 
the phase separation transition in a combined Falicov- 
Kimball-static-Holstein model. We see that increasing 
the electron-phonon coupling forces the system to be 
more homogeneous and suppresses the phase separation. 
We conjecture that a finite value of g will be sufficient 
to remove phase separation entirely, but are unable to 
access that region of the phase diagram due to numerical 
difficulties. Our results could have an impact on CMR 
materials that have both disorder and electron-phonon 
interactions. Any phase separation induced by the disor- 
der may be suppressed by the electron-phonon coupling. 
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